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We present a formulation of domain wall fermions in the Schroedinger functional by following the 
universality argument given by Liischer To check whether the formulation works, we examine the 
lowest eigenmode of the free domain wall fermion operator. We confirm that the theory belongs 
to a correct universality class and that the eigenvector is localized near the boundaries of the fifth 
dimension. We also investigate the chiral symmetry breaking structure of the four dimentional 
effective operator We observe that the bulk chiral symmetry breaking disappears for a large fifth 
dimensional size, while the breaking originated by the boundary effects persists and exponetially 
decays away from the time boundaries. 
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1. Introduction 

In the study of CP violation by CKM unitary triangle analysis, hadron matrix elements of four- 
fermion operators, such as Bk, play a vital role. Accurate calculations of this quantity from first 
principles are an important task for the lattice QCD community. In such calculations, having chiral 
symmetry is crucial to avoid an operator mixing problem which causes uncontrollable systematic 
errors. Therefore, the RBC/UKQCD collaboration [jl|] is currently using lattice chiral fermions, 
such as, domain wall fermions (DWF), to compute Bk- In the course of the computation, there are 
many sources of systematic errors which one has to control. Among them, the non-perturbative 
renormalization (NPR) could be a serious source. At the moment, the collaboration has been using 
conventional schemes, namely, the RI/MOM scheme and its variant [§, However, this scheme 
could be subtle depending on the target quantity. To avoid such difficulties, a new scheme was 
invented, known as the Schrodinger functional (SF) scheme [Q]. This scheme provides a reliable 
way of estimating errors in the NPR. If one wants to use this scheme for the renormalization of Bk 
given by the RBC collaboration, first of all, one has to formulate DWF in the SF setup. This is a 
purpose of this paper. 

In the SF setup, a formulation of lattice chiral fermions is in fact a non-trivial task. This can 
be seen as follows. First, let us consider the continuum theory with the SF boundary conditions, 

P+yif{x)l,=o = P-Y{x)\xo=t = 0, (1.1) 
V{x)P-\xo=Q = Vix)P+\xo=T = 0, (1.2) 

with P± = (1 ±Yo)/2. These boundary conditions break chiral symmetry explicitly. This fact 
indicates that the fermion propagator in such a theory does not anti-commute with 75 and breaking 
term is given by 

Y5D-'+D-'y5=A. (1.3) 

The breaking term A is supported at the time boundaries, whereas the anti-commutation relation 
holds in a bulk. This is a situation of the continuum theory. On the other hand, if one naively defines 
lattice chiral fermions by using the SF Wilson kernel operator, such a chiral operator automatically 
satisfies the Ginsparg-Wilson (GW) relation [0] and its inverse cannot reproduce the continuum 



results in eq.(1.3). This shows a difficulty of formulating the chiral fermions on the SF. 

To overcome such a situation, Taniguchi made the first attempt of formulating overlap 
fermions in the SF by using an orbifolding technique. Furthermore, he extended his idea to DWF 
[^] and then he and his collaborators obtained the renormalized Bk in the quenched QCD |^]. 
However, in dynamical simulations, his approach for overlap fermions has a sign problem for odd 
flavors. Although this problem can be removed for even flavors, the flavor symmetry is broken 



explicitly as an expense. Sint [ ]10| ] proposed a solution for even numbers of flavors which preserves 
the flavor symmetry exactly. However, this formulation can only be applied for the even number of 
flavors. This conflicts with the current trend toward large scale dynamical 2+1 flavor simulations. 
To overcome such a circumstance, Liischer [||] proposed a completely different approach following 
the universality and symmetry considerations. In that paper, he gave a formulation for overlap 
fermions. However, DWF from this approach has not been discussed so far, therefore we address 
this issue in this paper. 
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In the following, after introducing Liischer's argument, we give a formulation of DWF in 
the SF setup. And then we show some investigations to check how it really works. To enhance 
readability, in the following, the lattice spacing is suppressed unless necessary. 



2. Universality argument 

The argument in the continuum theory given in the previous section suggests that on the lattice, 
the GW relation has to be modified to correctly reproduce the continuum results. To this end, a chi- 
ral fermion operator also has to be modified compared with that of the usual infinite lattice. It is not 
so hard to break the GW relation itself and one can easily imagine that such a modification would 
be local and related with the boundary conditions. However, it is not clear how the SF boundary 
conditions come out? Liischer gave an answer to the question [^, which uses a symmetry consider- 
ation and an order counting. What he concludes is that the SF boundary conditions are natural and 
more stable than other boundary conditions in the continuum limit, because the dimension of the 
Dirichlet type boundary conditions is lower than that of the Neumann type boundary conditions. 
Bottom line is that it is guaranteed that the SF boundary conditions are automatically reproduced 
in the continuum limit without fine tuning if the lattice fermion operator is modified such that the 
GW relation is broken near the time boundaries. 

Before going into the concrete discussion of DWF, let us summarize a general instruction to 
formulate chiral fermions in the SF from the Liischer's argument. Firstly, including additional 
terms, which break the GW relation, to a lattice fermion operator such that the GW breaking is 
supported near the time boundaries. Secondly, keeping other symmetries, like C, P, T and flavor 
symmetries etc., which we do not want to lose. This instruction does not dictate a unique form of 
the operator, and actually there are an infinite number of possibilities. The important thing here is 
that the details of the operator actually do not matter as long as it fulfills these conditions. When 
these conditions hold, the corresponding lattice fermions are automatically guaranteed to follow the 
SF boundary condition in the continuum limit thanks to the order counting argument. By following 
this instruction, we define DWF in the SF in the next section. 



3. Universality formulation of DWF 

A possible form of the massless DWF operator in the SF is given by 

Dw + 1 -Pl cB 
-Pr Dw + 1 -Pl cB 
-Pr Dw + 1 -Pl + cB 

-Pr-cB D^ + l -Pl 
-cB -Pr Dw + 1 -Pl 

-cB -Pr Dw + 1 

This is a five dimensional block form and the block elements are four dimensional operator. 
Dwi^ns) is the Wilson fermion operator in the SF [|ri|] with the parameter nij and is supported 
in a range of the time 1 < ;co < T — 1. Here, the size for the fifth direction^ is taken as = 6 as 



D 



DWF ■ 



(3.1) 



We assume that L, is an even number. 
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an example. The additional terms which are proportional to the coefficient c are distributed along 
the cross diagonal elements. This is something like a generalized mass term and breaks the chiral 
symmetry explicitly. The block element B is now supported only near time boundary, 

B{x,y) = 5^,y5,,^y,Y5{5xo,iP- + Kt-iP+)- (3-2) 

This time dependence comes from a motivation that we do want to maintain the chiral symmetry 
in a bulk in time. We will check that the breaking of the GW relation is really supported near the 



time boundaries exponentially (see Section |4-.2| ). In this way, such fifth dimensional coordinate and 
time dependences are determined. Next, how about the spinor structure? Actually, this structure 
may be fixed to some extent by requiring the C, P, T symmetries and the Fs-Hermiticity. These 
requirements are not so strong to determine the spinor structure completely and there is some 



freedom. The structure given in eq.(|3.2D is one of many solutions. Actually, we examined several 
choices of the spinor structure in the boundary term and confirmed the universal results in the 
continuum limit for the lowest eigenvalue. 

The coefficient^ c 7^ is a parameter which can be used to achieve the 0(a) improvement. The 
physical mass term can be added in the same way as the infinite lattice case. 

4. Check of formulation 

The next step is to check how it works. We investigate several things as follows. First, we study 
the lowest eigenmode for the free operator to see the importance of the presence of the boundary 
term. Second, we address the breaking of the GW relation for the four dimensional effective 
operator. This is to see what the chiral breaking term looks like. We also checked scaling properties 
of the free spectrum of the lowest ten eigenvalues and carried out the one-loop calculation and 
obtained sensible results. However we do not show them here due to limited space. In the near 
future, we will show them in a full paper. 

4.1 The lowest eigenmode for free massless operator 

Let us see the lowest eigenmode of the free massless operator D^^pDowF 

^DWF^DWF Y = KimW- (4- 1 ) 

To show the importance of the boundary term, we consider two cases with c = 1 and c = 0. In 
the free case, the Fourier transformation can be performed for the spatial directions and we project 
onto the zero momentum configuration. 

For c = 1, we obtained the minimum eigenvalue, T^Xmin = 2.59... This value is rather close to 
the continuum one 71^/4 = 2.47.. in Ref.[12]. The difference is considered to be a lattice artifact. 



The eigenfunction in Figure |I] (left panel) shows the expected behavior, namely, being localized 
in the fifth direction while propagating in the time direction. This is apparently a physical mode. 
On the other hand, when switching off the boundary term c = 0, the eigenvalue 0.61... is quite far 
from the continuum value. Actually, a scaling study with larger lattice sizes shows that the theory 



: is not allowed as we will see in Section |4. l| . 
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Figure 1: ||v^(.5[:o,s)||spin with the zero spatial momentum with the parameters T = Lg = 40, = and 
m^ — l. The norm for the eigenvector is taken for the spinor space. The left (right) panel is for c = 1 (c = 0). 



does not belong to the correct universality class. Furthermore, the eigenfunction in the right panel 
of Figure [T] appears to be propagate in the fifth direction. This is not a physical mode anymore. 

This result shows that the presence of the boundary term is essential to correctly produce the 
continuum results. In other words, this term is essential to be in a correct universality class. 

4.2 Chiral symmetry breaking by boundary effect 

Next, let us check that DWF in the SF has an expected chiral symmetry breaking structure. To 



this end, for the effective four dimensional operator 1 13, 14], 

det[D^^')] = det[DDWF/£»Pv] , (4.2) 
let us consider the breaking of the GW relation 

A(^') = 75d(^') +Di^'V5 -2d1^:V5£»S'^ (4.3) 

For finite L^, there are two sources of chiral symmetry breaking, namely, a bulk source and a 
boundary source. The former can be removed by taking to infinity. After this limit, the remaining 
chiral symmetry breaking must be due to the boundary source only. Actually, it is known that the 
corresponding breaking for overlap fermions exponentially decays away from the time boundaries 
[^. We want to check this kind of phenomena for DWF in a free case numerically. 

Figure ^ shows the magnitude of the chiral symmetry breaking A'^^^^ in eq.(43) with the zero 
spatial momentum configuration. The norm in ||A||spin is taken for the spinor index. For = 4, 
we can see large chiral symmetry breaking in diagonal elements not only near boundaries but also 
in a bulk. However, at increased L^, the bulk breaking gets smaller and then finally it disappears 
for Lg = 24. In the end, the remaining breaking is exponentially localized near boundaries and we 
observe the expected behavior mentioned above. 
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Figure 2: log | |A*^^'^ (-^cyo) I Ispin for the zero spatial momentum with the parameters T = 20, mf = 0, = 0, 
= 1 and c = 1 . 



5. Concluding remarks and outlook 

By following the universality argument, we constructed DWF in the SF setup. And then we 
checked some fundamental properties of the theory, namely, that the lowest eigenmode is local- 
ized around the fifth dimensional boundaries and the GW relation breaking for the effective four 
dimensional operator has the expected structure. Although we do not show it in this paper, we also 
checked the universahty at the one-loop level. In total, our formulation works very well. 

We are now ready to study NPR of Bk for any flavors with DWF by using the SF scheme. 
Finally, we comment that by using our formulation, one can carry out a search for a conformal 
window. Since chiral fermions have benefits such as requiring no mass tuning to realize the mass 
independent renormalization scheme and no constraint on the number of flavors, our formulation 
can play a crucial role for solid quantitative studies. 

S.T would like to thank Sinya Aoki, Yasumichi Aoki, Norman Christ, Michael Endres, Taku 
Izubuchi, Changhoan Kim and the members of the RBC collaboration for helpful discussions. This 
work is supported by the U.S Department of Energy under Grant No. DE-FG02-92ER40699. 



References 

[1] RBC-UKQCD Collaboration, C. Allton et. al. Physical Results from 2+1 Flavor Domain Wall QCD 



and SU(2) Chiral Perturbation Theory, Phys. Rev. D78 (2008) 1 14509 [ |0804 . 0473[ |. 
[2] Y. Aoki et. al., Non-perturbative renormalization of quark bilinear operators and Bk using domain 



wall fermions, Phys. Rev D78 (2008) 054510 [ |0712 . 106l[ ]. 
[3] C. Sturm et. al., Renormalization of quark bilinear operators in a MOM- scheme with a 



non-exceptional subtraction point, Phys. Rev. D80 (2009) 014501 [0901 . 2599| 



6 



A formulation of domain wall fermions in the Schrodinger functional 



Shinji Takeda 



[4] M. Liischer, R. Narayanan, P. Weisz and U. Wolff, The Schroedinger functional: A renormalizable 
probe for nonabelian gauge theories, Nucl. Phys. B384 (1992) 168-228 [ |hep-lat/92 07 9[ ]. 

[5] M. Liischer, The Schroedinger functional in lattice QCD with exact chiral symmetry, JHEP 05 (2006) 

042 [ ]hep-lat/0603029[ ]. 

[6] P. H. Ginsparg and K. G. Wilson, A remnant of chiral symmetry on the lattice, Phys. Rev. D25 (1982) 
2649. 

[7] Y. Taniguchi, Schroedinger functional formalism with Ginsparg-Wilson fermion, JHEP 12 (2005) 037 

[ [hep-lat/0412024| ]. 

[8] Y. Taniguchi, Schroedinger functional formalism with domain-wall fermion, JHEP 10 (2006) 027 

[ |hep-lat/0 604 002| ]. 

[9] CP-PACS Collaboration, Y. Nakamura, S. Aoki, Y. Taniguchi and T. Yoshie, Precise determination of 
Bk and right quark masses in quenched domain-wall QCD, Phys. Rev. D78 (2008) 034502 
[ |0803.2569| ]. 

[10] S. Sint, The Schroedinger functional with chirally rotated boundary conditions, PoS LAT2005 (2006) 

235 phep-lat/0511034p . 

[11] S. Sint, On the Schroedinger functional in QCD, Nucl. Phys. B421 (1994) 135-158 
[ [hep-lat/9312079l ]. 

[12] S. Sint and R. Sommer, The running coupling from the QCD Schroedinger functional: A one loop 
analysis, Nucl. Phys. B465 (1996) 71-98 [ |hep-lat/9508012[ ]. 

[13] Y. Kikukawa and T. Noguchi, Low energy effective action of domain-wall fermion and the 



Ginsparg-Wilson relation, hep- 1 at/ 9 902 022 



[14] Y. Kikukawa, Locality bound for effective four- dimensional action of domain-wall fermion, Nucl. 
Phys. B584 (2000) 51 1-527 [|hep-lat/'9912056n. 



7 



